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Abstract 

- A theoretical and experimental study was undertaken into the crushing 
behavior of axially compressed short thin-walled open-section columns. The 
effect of initial geometry of panels as well as distribution and magnitude of 
shape imperfections on the efficiency of energy absorption was examined. 

Results of model tests on 0.1mm thick aluminum foil specimens have shown that 
the panels collapsing in the symmetric and asymmetric deformation mode provide 
respectively, upper and lower bound for the energy absorbed in any other buck¬ 
ling mode. In both of the extreme cases, the crush response of the panel was 
predicted theoretically with a reasonable accuracy. It is shown that an opti¬ 
mum design of columns against crush can be achieved by introducing a beneficial 
geometric imperfections of a specified magnitude so that the structure will be 
forced to collapse in the most energy efficient deformation mode. 
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Introduction 

The problem of reducing the weight of a vehicle body without 
sacrificing its strength, stiffness and crashworthiness properties 
has become a great concern to auto manufacturers. The optimization 
of thin-walled structures for energy absorption can be approached 
in many different ways. Weight reduction can be achieved by con¬ 
sidering alternative materials such as aluminum, high strength steel 
or fiber reinforced composites,[1],[2]. Another possibility is to 
strengthen properly the shell by means of stiffeners in order to 
impose a desirable and energy efficient deformation mode. One 
could also incorporate into the design initial imperfections of a 
specified magnitude which again are aimed at triggering the right 
buckling mode and forcing the structure to collapse in that mode. 
Finally, great potential exists in properly shaping the sheet metal 
structures for maximum crash resistance under prescribed geometrical 
constraints. 

The problem of shape optimization of thin-walled structures 
subjected to quasi-static crushing can be formulated as follows: 


For a given weight and prescribed total crush 
distance,maximize the energy absorption of the 
structure. 
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Dual formulation: 

For a given amount of the energy to be absorbed 
and for a prescribed total crushing distance, 
minimize the weight of the structure. 

Clearly, the application of the optimal criterion requires 
a careful determination of the admissible class of structures and 
boundary conditions. Very little work on the more fundamental nat¬ 
ure has been done in this new area. It is generally understood that 
closed - section members under compressive forces are far more 
efficient from the point of view of energy absorption than o^en - 
section structures. It has also been observed that the crushing 
force increases with the number of sharp corners (angles) of pris¬ 
matic or channel section tubes [3]. According to various authors, 
the mean crushing force is proportional to the gauge thickness h 
raised to the power 1.5 f 2. On the other hand, only slight de¬ 
pendence of the force was observed in the circumference of the 
tube. These findings gave rise to the concept of using long, pro¬ 
gressively collapsing tubular members for crash protection of auto¬ 
mobile bodies. Its limitation is set by the condition of overall 
loss of stability which restricts the length to width ratio Vb 
to relatively low values. On the other hand, closed - section mem¬ 
bers with small Ji/b have to be excluded from most of the designs 
because of package requirements. In order to get the required 
crash resistance, one has thus to turn to open section shells. 

In the present paper, the response of such shells under quasi¬ 
static compressive loading is studied to get some insight into the 
extremely difficult problems of the formation and motion of the 
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plastic folds and wrinkles. Rather than solving the optimization 
problem for general shells, we shall restrict the analysis to 
axially loaded open panels forming one quarter of prismatic thin- 
walled tubes. The dependence of the crushing characteristics of 
the panels on imperfections, boundary conditions and the transition 
radius was studied experimentally and theoretically and important 
parameters were identified which are controlling the process of 
energy absorption. 

Experimental Procedure 

Fourty-three model tests were run in three series on 0.1 mm 
thick aluminum foil specimens with dimensions shown in Figure 1. 

The height of all specimens was the same, Z = 75 nun; the ratio 
Z/h being typical to large panels of the automobile body such as 
the hood, fenders, trunk lid, etc. The total circumference of the 
specimens was also held constant c = 100 mm but the shape varied from 
the symmetric angle element to the cylindrical panel with the cen¬ 
tral angle equal to The considered transition radii were respec¬ 
tively, p = 0,5,20,25,35 and 65 ram. The lateral surface of the 
specimens was fixed to the bottom and top plates by means of a 
Scotch tape whose bending stiffeness was much smaller than that of 
the foil. Thus, the boundary conditions were of the simple sup¬ 
ported type along the horizontal edges and free along the vertical 
edges. Intentionally, little care was given to eliminate the shape 
imperfection in order to best resemble the real world situation. 

On a certain number of the specimens, intial imperfections were im¬ 
posed described by the amplitude £, half wavelength X and number 
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of waves n. Figure 2. All specimens were crushed up to one half 
of its original length in a kinematically driven testing machine. 
Both the bottom and top plates were fixed to the cross-head plates 
to prevent rotation. In each experiment, the force-shortening 
characteristics was recorded on the x-y plotter. From these dia¬ 
grams, the energy absorbed at 30 mm crush and the average force 
level were calculated. Photographs of all specimens were taken 
before and after each test. In the first series of tests, "perfect" 
structures were crushed with the transition radius p as a design 
variable. In the second series, the effect of large initial imper¬ 
fections was studied for specimens with selected transition radii. 
Finally, the effect of the deformation mode on the energy absorption 
of panels with p = 25 nun was examined in the third series. 

Experimental Results and Discussion 

A typical load-deflection relationship is shown in Figure 3. 

One or several peaks of the load are observed, the first being usu¬ 
ally, but not always, the maximum one. Each distinct peak is assoc¬ 
iated with the formation of a different local buckling mode. The 
force is generally diminishing with deflections. In most cases, the 
pattern of folds and wrinkles is very complicated and depends on 
initial random imperfections. In spite of this seemingly unpredic¬ 
table behavior of panels, some interesting features of the crush¬ 
ing process have been identified. 

Effect of the Panel Curvature 

The total energy absorbed at 30 mm. crush distance was plotted 
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against the transition radius p for all test pieces without imposed 
imperfections (Figure 4). While the scatter of experimental points 
for each value of the radius p is quite significant, a well defined 
tendency can be observed. The full line in Figure 4 represents the 
best fit of all experimental points? Thus, an optimum structure con¬ 
sists of a cylindrical panel supported on both sides by two flat 
plates. The energy absorption of the panels with the largest radius 
p=65mm was only one half of the optimum value while the efficiency 
of angle elements with p = 0 falls between these two extreme cases. 

Effect of Imperfections 

The objective of the crushing tests on imperfect structures 
with imposed large initial imperfections was to see whether it is 
possible to excite deformation mode with short wavelength and get 
a progressive collapse of the considered open-structure members. 

In most cases we have introduced very large imperfections, sometimes 
exceeding by two orders of magnitudes, the gauge thickness of an 
element! Panels with large imperfections generally followed the 
collapse mode imposed by imperfections regardless of their wave¬ 
length. In constrast to the closed-section members, the deforma¬ 
tion mechanism was not progressive; we have simply observed a sim¬ 
ultaneous folding of the material along the stationary hinge line. 

The associated force - deflection relationship was relatively smooth 
but far below the corresponding curve for "perfect" shells. The 
dramatic decrease in the maximum force is mainly attributed to the 
large magnitude of initial imperfections. Examples of crushed spe¬ 
cimens with large symmetric and asymmetric imperfections are shown 

*It has clearly a maximum for p = 25 mm 
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respectively, in Figures 5 and 6. Specimens with much smaller ini¬ 
tial imperfections behave in a different manner, as exemplified in 
Figure 7. When the wavelength of imperfection is relatively short, 
the structure disregards its presence and eventually assumes an 

l 

asymmetric rather than symmetric deformation mode with H = The 
associated reduction in the force level is now much less pronounced. 
Some aspects of the transition from the symmetric to asymmetric buck¬ 
ling mode will be discussed further in the paper. 

The third series of tests was run on panels with an optimum 
radius, p = 25mm and no imposed initial imperfections. The random im¬ 
perfections cause the structure to assume in each case a different, 
and in most cases, quite a complicated collapse mode leading to dif¬ 
ferent load-deflection characteristics. However, two type of acti¬ 
vated modes are of particular interest. One is a regular asymmetric 
o 

mode with H = j. Figure 8, while the other one is a symmetric mode 

£ 

with H = Figure 9. The plastic energy dissipated in these modes 
were found to furnish respectively, the lower and upper bound on the 
energies dissipated in all other modes . This is the most important 
conclusion of the present experimental program. These findings bear 
important implications for the analysis and design of the thin pan¬ 
els. Now, the worse and best possible case of energy absorption can 
be predicted on the basis of a relatively simple collapse mechanism. 
Consequently, the analyst should not be worried about the effect of 
random imperfections and about the difficulties associated with con¬ 
sidering in the calculations more complicated buckling pattern. On 
the design side, the problem of maximizing the energy absorption can 
be reduced to the problem of triggering the desirable symmetric 











buckling mode. 

In the remainder of the paper, attention will be focused on 
the plastic crushing analysis and elastic post-buckling analysis 
of panels with symmetric and asymmetric collapse mechanism. 


General Analysis of the Plastic Crushing Process 

Consider a typical fold line which moves down the shell ele¬ 
ment as the deformation progresses goes on, Figure 10. Since the 
fold line always lies in a certain plane, (see Ref. [4], it is reason¬ 
able to model the deformed shell element by a section of an elliptic 
toroidal surface, four sections of a cylindrical surface and four plane 
trapezoidal elements. Such a model, first suggested in [5] and [ 6 ], 
is fully consistent geometrically and kinematically. In previous 
models considered in the literature, the toroidal surface was not 
considered which leads to the discontinuities in the displacement 
field at the corner point C, [7], [ 8 ], [9]. To simplify the calcula¬ 
tions, the elliptical toroidal surface will be approximated by a cir¬ 
cular one defined by a central angle 2 g and a small and large radii, 
denoted respectively by r and R, Figure 12. The global geometry 
of the crushing process is described in Figure 11. The angle 
between two adjacent sides and the total length of segments AC+CD 
are kept constant while the angles y (a) andTr-g(a) diminish with the 
parameter of the process a. 

The angles a, 6 , y and 41 are related by (tg stands for tangent). 


tg Y 


-tai 

sing 


2 

cos 8 


1 -sin a 


1 + ( 


sina 

tg<J> 


2 

) 


( 1 ) 


The crushing distance 6 and crushing velocity are related to the 
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wavelength H and a through 

<5 = 2H (1-cosa), 6 = 2H(sina)a 


( 2 ) 


According to the assumed model, four mechanisms contribute to 
the energy dissipation. There are: 

1. Extension in the hoop direction of the toroidal surface. 

2. Change of curvature in the circumferential direction from 
positive 1/R to negative - 1/R. 

3. Bending and rebending along moving plastic hinges 1 and 2. 

4. Pure bending in stationary plastic hinges 3 and 4. 

The strains developed in the toroidal surface may sometimes be 
very large thus requiring consideration of the variable thickness. 
The exact expressions for the dissipated energies in each of these 
mechanisms was calculated in {53 and [6]. In the present analysis, 
we shall use the following simplified formulas, valid for constant 
thickness: 
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component of the velocity of the line AC or CD is 


V - HCcoscOa , v t - ^ 


Equating the sum of (3) 7 ( 6 ) to the rate of work of external 
forces 

*ext - Pi <8) 

one gets a final expression for the magnitude of an instantaneous 
axial force necessary to maintain the plastic flow, 


M o = se h coa * ctga + 2a E ^r + 


H ctga 1 + 2 — ~b _ 

r siny tgip H sing 


We shall consider now two special cases of Equation 


Crush Prediction in the Asymmetric Mode 

Consider a panel with p = 0, in which an asymmetric mode is 
activated. The angle 4) is equal to n /^r so that the formula (9) 


is reduced to 


= 5.66 ^ ctgg + 2 ctgg + — fctgg + 2 pj- 


C 1 


where 


f(a) = • 
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The right hand side of (10) can be formally minimized with respect 

3 P 

to the radius r and wavelength H. The conditions /. =0 and 

o X 

3p / 3h = 0 yield 


r 

opt 


0.42 



(ID 


H o P t 


0.89 



_C 

/* f 


IT 
cos a 


( 12 ) 


Thus, the optimum values of r and H turn to be a function of a 
which is inconsistent with the assumed mechanism of deformation 
in which r and H were both taken as constant. While the variation 
of r with a can easily be accommodated by the present solution with¬ 
out requiring much of the additional energy to be dissipated, the 

8 P 

consequences of the condition, / 3H = 0 should be discussed 
in more detail. 

— 31 5 

The variation of dimensionless wavelength, H = H/\JC h with 
the parameter of the process shown in Figure 13. The wavelength 
of the corresponding elastic buckling mode is H e ^ = -j-, which 
yields H = 3.75 (broken line in Figure 13). Experiments show 
that the panels collapse in the asymmetric mode with the elastic 
wavelength. A considerable variation of H, predicted by (12) would 
not be possible without dissipating additional energy on the down¬ 
ward motion of the otherwise stationary hinge lines ACD and corner 
point . The necessary calculations have not been performed but it 
is anticipated that the new energy terms would prevent any major 
variation of H. It is thus, very likely that the wavelength of 
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the plastic collpase mode is decided early in the deformation pro¬ 
cess and a possible adjustment of the H takes place in the stage 
of an elastic post-buckling behaviour of the panel. 

Substituting (11) and H = 37.5 mm, c = 100 mm, h = 0.1 mm 
into (10), we get the parametric representation of the force-dis¬ 
placement relationship of the considered panel 


— = 192 + 2a^ ctga + 5.3/sina (13) 

6 = 2H -^l-cosa^ 

The prediction of (13) has been compared with P - 6 curve measured 

in tests No. 32 and 35. The agreement is very good except for the 

initial phase when the present theory is not valid. (Figures 14 and 

15). In matching the theory with experiments, we have assumed the 

KG 

yield stress of the aluminum foil to be a =8 / 2 . 

•* o mm 


Crush Prediction in the Symmetric Mode 

The simplified model of the symmetric collapse mode is shown 
in Figure 16. The existence of such a mode can be easily evidenced 
by making an inextensible paper model of the panel. Now, there 
are two hinge lines and hence, all terms in the energy balance 
equation, except the contribution of the stationary hinge lines, 
should be doubled. The angle tp becomes Vg and from (9) one can 
get an expression for the crushing force 


pp = 14.78 yj ctga + 9.65actga + 4.82^fctga + 3.14 


(14) 
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A minimization of the right hand side of (14) with respect to r 
and H gives 


r opt = °' 57 -\1 


hHf 0 


-1 


(15) 


H , = 0.52 
opt 


A 


fBcos a 


(16) 


The plot of Eq. (16) reveals that the predicted wavelength of the 

symmetric mode is two times smaller than that corresponding to the 

asymmetric mode for all values of a, Figure 13. This tendency is 

o 

confirmed by test No. 30 where H was of the order of / rather than 
£ 

2 . Substituting (15) back into (14) , we obtain the sought for 
expression for the crush resistance of the panel 


P 

M 

o 


16.85. 


H 




j fM 9.65a^ ctga + 


3.14 


C 1 
H sina 


(17) 


Plots of Eq. (17) for several chosen values of H are shown in 

Figure 17. The correlation with the experimentally measured force 

is not as good as in the case of the asymmetric mode of deformation. 

The present solution correctly predicts the general character of 

the P - <S relationship but overestimates the force level for a 

£ 

realistic value of the wavelength H = j = 18.75. It should be 
noted that the activation of the symmetric mode, clearly visible 
in Figure 6 is preceded by the development of some other modes 
near the bottom plate, which explains a poor agreement of the 
theory with test results at the initial stage of the crushing pro- 
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cess up to 6 = 10 mm. Some discrepenci.es between the theory and 
experiments may also be attributed to the fact that the calcula¬ 
tions were made for panels with p = 0 while the measurements were 
taken for p = 25 mm. 

In the same figure shown is also the P-6 curve corresponding 
to the asymmetric mode (broken line). Under any circumstances, 
this curve is much lower than those representing a symmetric mode 
solutions and so is the dissipated energy. 

In summing up, it can be stated that the present solution 
gives a good prediction of the crushing resistance of panels de¬ 
formed in an asymmetric and symmetric mode, provided that the wave¬ 
length of the collapse mode is properly determined. The energy 

absorption calculated in these two modes furnishes 
a lower and upper bound on energies associated with any other de¬ 
formation modes of panels with various p and random imperfections. 



(18) 






where D is the bending rigidity of a plate, b - its width and 
the coefficient k depends on the aspect ration b/a which in our 
case is With the present dimensions of the plate c = 100 mm, 

£ = 75 mm. This coefficient is equal to 


[0.9 


k =< 


asymmetric mode 


1.34 


symmetric mode 


The buckling mode corresponding to the lowest buckling force is in 
both cases a half of the sinusoidal wave, i.e. H sfA • After the 
critical load is reached, a post buckling path is followed 
for both symmetric and asymmetric case. Figure 19. Thus, a per¬ 
fect elastic angle element will always buckle in an asymmetric mode 
which requires a much lower critical load. 

The situation might be different for an imperfect structure 

since now the two modes may interact with each other. In the 

majority of cases, the most degrading effect on the load have the 

imperfection in the shape of the buckling mode [11]. Denote by y 
___ a 

and ^ , the maximum amplitude of initial imperfections respec- 

s __ __ 

tively in the asymmetric and symmetric mode. If ^ ^ 0 an< * t = 0, 

a s 

the structure will follow the equilibrium paths in the asymmetric 
mode, Figure 19. If the compression of the plate is continued, 
this will eventually lead to the lowest possible energy dissipa¬ 
tion. Suppose now that f = 0 and £ f 0. The structure will ini- 
tially follow the symmetric mode deformation but at a certain point 
such as F may bifurcate into the asymmetric mode. The question 
which remains to be answered is whether this will happen and if so, 
when and how the coordinates of the point F depend on the magnitude 
of the initial imperfections t g . 
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The important problem is being studied using Koiter's theory 
of post-buckling behaviour of imperfect structures with mode inter¬ 
actions, [12],[13] and the results will soon be published, [14]. 

If we succeed in controlling the deformation process so that the 
structure will not pop back from the symmetric mode to the pref¬ 
erable (giving lower force) asymmetric one, then the way for opti¬ 
mizing the structure for crash would be open. 

Concluding Remarks 

We have shown that the energy absorption of compressed panels 
can be considerably increased if a proper shape of the cross-section 
is assumed. Within the considered class of structures, the best 
shape consists of a cylindrical panel representing one quarter of a 
complete cylinder supported on both sides by flat plates. The 
presence of random initial imperfections give rise to quite a 
variety of collapse mechanisms. However, two specific types of im¬ 
perfections are of utmost importance. It was found that the imper¬ 
fections in the asymmetric and symmetric elastic buckling mode have 
respectively the most degrading and beneficial effect on the energy 
absorbed by the panel. In both of the extreme cases, the crash 
response of the panel can be predicted with a reasonable accuracy. 
Investigation is under way regarding the stability of the buckling 
modes imposed by symmetric imperfections. A new series of experi¬ 
ments is also being planned to demonstrate the effect of the bene- 
fitial imperfections on the crush response of the panels. 

The required symmetric collapse mechanisms may be triggered 
in many different ways, for example, by attaching string stiffeners 
with initial curvature to the free edge of the panel. The problem 









- 16 - 


of interaction of stiffeners and a shell with the view of maximizing 
the energy absorption will be the subject of a future research. 

It is believed that some of the factors identified in the pre¬ 
sent study would also affect the energy absorption of sheet metal 
structures with different geometry and boundary conditions. 
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FIGURE 19 
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